Stationary state volume fluctuations in a granular medium 
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A statistical description of static granular material requires ergodic sampling of the phase space 
spanned by the different configurations of the particles. We periodically fluidize a column of glass 
beads and find that the sequence of volume fractions (j> of post-fluidized states is history indepen- 
dent and Gaussian distributed about a stationary state. The standard deviation of cj> exhibits, 
as a function of (f>, a minimum corresponding to a maximum in the number of statistically inde- 
pendent regions. Measurements of the fluctuations enable us to determine the compactivity X, 
a temperature-like state variable introduced in the statistical theory of Edwards and Oakeshott 
[Physica A 157, 1080 (1989)]. 

PACS numbers: 45.70.-n, 05.40.-a, 64.30.+t, 47.55.Kf 



Granular materials consist of a large number N (typi- 
cally more than 10 6 ) dissipative particles that are massive 
enough so that their potential energy is orders of magni- 
tude larger than their thermal energy. The large number 
suggests that a statistical description might be feasible. 
Edwards and coworkers developed such a description 
with the volume V of the system, rather than the en- 
ergy, as the key extensive quantity in a static granular 
system. The corresponding configuration space contains 
all possible mechanically stable arrangements of grains. 

Brownian motion is insufficient for a granular system 
to explore its configuration space, so energy must be sup- 
plied by external forcing such as tapping [2j , shearing |3j 
or both The theory of Edwards requires that the 
forcing assures ergodicity: all mechanically stable config- 
urations must be equally probable and accessible. A nec- 
essary condition for ergodicity is history independence: 
the physical properties of the system must not depend 
on the way a specific state was reached. History in- 
dependence has previously been demonstrated only by 
Nowak et al. for tapped glass beads at volume frac- 
tions 4> > 0.625. 

In this paper we explore the configuration space using 
a periodic train of flow pulses in a fluidized bed. A sta- 
tionary column of glass beads in water is expanded by an 
upward stream of water until it reaches a homogeneously 
fluidized state [5j , and then the flow is switched off. The 
fluidized bed collapses and forms a sediment of volume 
fraction <f>, which we find depends in a reproducible way 
on the flow rate Q of the flow pulse. This forcing results 
in a history independent steady state where the volume 
exhibits Gaussian fluctuations around its average value. 

A central postulate of the Edwards theory is the ex- 
istence of a temperature-like state variable called com- 
pactivity X — dV/dS. The entropy S is defined in anal- 
ogy to classical statistical mechanics as S(V, N) = X In Q, 
where is the number of mechanically stable configura- 
tions of N particles in V, and A is an unknown ana- 
log to the Boltzmann constant. The assumption that 



X is a relevant control parameter in granular systems 
has found support in simulations of segregation in bi- 
nary mixtures J3 , compaction under vertical tapping || , 
and shearing '9]. However, no measurements of X have 
been reported. In this paper we determine X from the 
measured volume fluctuations using a method suggested 
in 0. 

Experiment. The apparatus is shown in Fig. ^a). 
A square bore glass tube (24.1mm x 24.1mm) con- 
tains about 3.6 x 10 6 beads (soda-lime glass, d — 250 ± 
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FIG. 1: (a) Experimental setup. The inset shows the cor- 
relation between the height changes Ah during a flow pulse 
to the fluidized bed, as seen by the two cameras (800 flow 
pulses, Q — 60mL/min). (b) Development of the bed height 
during flow pulses of Q = 40 mL/min. The arrows (f ) indicate 
measurements of h. 
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FIG. 2: History independent steady state: (a) The volume 
fraction <f> of the sedimented bed reaches a steady state that 
depends only on the flow rate Q. (b) The inset shows the in- 
crease in volume fraction as Q was decreased in 14 steps from 
60mL/min to 23mL/min, followed by a decrease in volume 
fraction as Q was increased again; each step consisted of 50 
flow pulses. The main graph (b) shows that the steady state 
volume fraction is the same for increasing (o) and decreasing 
(•) Q; the dashed line is a fit to |T). 



13 /xm, density 2.46g/cm 3 , MO-SCI Corporation). The 
beads are fluidized with pulses of temperature-controlled 
(23.0 ± 0.1 °C) de-ionized water. 

Flow pulses were generated by a computer-controlled 
syringe pump. During a flow pulse of length r pu i sc 
(Fig Gib)), the bed expanded until its height reached a 
stable value where the time-averaged viscous drag force 
on each particle corresponded to its weight; r pu i sc was 
chosen so <3t pu i sc = 30 mL. After each flow pulse, the 
bed was allowed to settle into a mechanically stable con- 
figuration, which was determined by measuring the time- 
resolved correlation of laser speckle patterns |l0| . A wait- 
ing time T wa it of 30 s was found to be ample for achieving 
a jammed state ^lj- After r wa it the volume fraction <f> 
was determined by measuring the bed height h with two 
CCD cameras at a 90° angle. Averaging over the width of 
the bed yielded height values with a standard deviation 
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FIG. 3: Volume fraction fluctuations around the steady states 
shown in Fig. |5Ja) are Gaussian distributed. The upper row 
shows examples of the volume fraction measurements as a 
function of the number of flow pulses; </> aV g is indicated above 
each column. The lower row shows the corresponding proba- 
bility distributions obtained from 800 flow pulses; the curves 
are Gaussian fits. 



of only 0.008d for a bed at rest (cf. inset of Fig.QJa)). 

History independent steady state. Starting with an ini- 
tial <f> — 0.581 prepared with a single flow pulse of Q 
= 60mL/min, the volume fraction quickly approaches a 
Q-dependent constant value </> avg , as shown in Fig [2] (a). 
Then with successive flow pulses, (f> fluctuates about (fravg- 
The history independence is demonstrated by ramping up 
and down in flow rate (Fig.^b)); <p depends only on Q 
of the last flow pulse, not the earlier history of the bed. 
The remaining differences in </> avg for increasing and de- 
creasing Q are on average only 1.6 x 10 ~ 4 . These small 
variations are correlated with the viscosity changes due 
to the temperature drift of ±0.1 °C during the course of 
the experiment (35 h). 

For slow sedimentation, corresponding to Q — + oo, <f> 
should converge to its random loose packing value, <^rlp 
[l2T |. In the absence of a theory, we fit the data to a 
phcnomenological equation, 



<KQ) = <t> 
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(Fig. □»). This yields </> RL p = 0.573±0.001 (o= 
0.365 mL/min, 6=12.9mL/min), where the uncertainty 
arises not from the very small uncertainty in bed height 
measurement (corresponding to 0.0003 in <fi) but from 
variations in the cross section of the square tube. 

Volume fluctuations. Gaussian fluctuations in </> are 
observed for the whole range of (f> studied; examples are 
shown in Fig. The standard deviation of these fluc- 
tuations rj0, like 4>avg, was found to be history indepen- 
dent. The variation of as with cj) fits a parabola with 
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FIG. 4: (a) Volume fraction fluctuation measurements fit a 
parabolic curve, (b) number of particles in a statistically in- 
dependent region, and (c) compactivity as a function of the 
average volume fraction. At ^ m the system has a maximum 
in the number of statistically independent regions (cf. ©). 



a minimum at <p m — 0.596 (Fig. E{a)). The minimum 
corresponds to a maximum in the number of statistically 
independent spatial regions Ni at the moment of solidifi- 
cation, 




where k w 0.2, estimated by Nowak et al. [2j from the 
maximal volume fraction fluctuations a single spatial re- 
gion could undergo. Using and the parabolic fit in 
Fig.^Ja), we estimate the average number of particles per 
statistically independent region, N ppr = N/Ni (where TV 
is the number of particles). Figure 0Jb) shows N ppr to 
be 40 for random loose packing, it reaches a minimum 
of about 1.8 at <j> m and then increases again. This con- 
trasts with the roughly constant value of N ppT found in 
previous experiments J2J for <f> in the range 0.625-0.634. 

We suggest that the minimum of iV ppr at cp m corre- 
sponds to a crossover of competing mechanisms. With 
increasing <f> above 4>RLP, where force chains stabilize 
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FIG. 5: Influence of particle surface roughness. On the left 
are scanning electron microscopy images of the surface of new 
and aged glass beads. The plot on the right shows the vol- 
ume fraction fluctuations of new beads (■) together with a 
parabolic fit (solid line); to minimize aging, measurements 
were made with only 200 flow pulses. The dashed line is the 
parabolic fit for aged beads from Fig. [I] 



the medium, the number of contacts per grain increases 
[l3|. and the probability that the displacement of one 
grain will mechanically destabilize other grains decreases. 
Consequently N ppl decreases. But with increasing <j>, the 
free volume per particle decreases. The latter effect is 
dominant for (f> > <f> m , where geometrical constraints al- 
low movement only as collective process of an increasing 
number of beads. Hence N ppl - and increase. At the 
maximum random jammed packing, 0mrj = 0.64, any 
motion would require a rearrangement of the whole sys- 
tem [y|. 

Knowing the dependency of on cf) enables us to de- 
termine Edward's compactivity X using a granular ver- 
sion of the fluctuation dissipation theorem derived by 
Nowak et al. Q ■ We obtain 

m Arlp X ( ( t ) ) 

where we have used -X"(<^rlp) = oo from granular statis- 
tical mechanics, and p is the particle density. We solve 
P)l numerically by substituting o$ from the parabolic fit 
(Fig. Ufa)) and starting at </>rxp determined from the fit 
with . This first experimental result obtained for X is 
shown in Fig. Ufc); the values involve the unknown mul- 
tiplicative factor A. No theoretical or numerical work has 
predicted how X depends on cp. Knowing X will allow 
further tests of granular statistical mechanics. 

Aging and friction. Our experiments revealed an ag- 
ing effect due to a change in the frictional properties of 
the beads. Friction between beads decreases with multi- 
ple collisions, which lead to a decrease in the number of 
surface asperities, as illustrated in Fig. The data we 
have presented show no aging during the 12,000 pulses 
in the experiment; the beads had been previously used 
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for more than 45,000 flow pulses. Just as for aged beads, 
the volume fraction standard deviation for new beads 
fits a parabola, as Fig. [S] illustrates; the minimum is at 
</> m =0.587. The shift of the minimum to a value lower 
than lower m =O.596 for aged beads is consistent with 
our argument about statistically independent regions: a 
higher values of /x makes dilute states less fragile and low- 
ers the onset of the necessity of collective motion. For 
new beads a fit to 0} results in </>rlp = 0.568 ± 0.001, in 
accord with </>rlp = 0.566 ± 0.004 interpolated from the 
measurements (presumably with new beads) in 1 1 i2t | . This 
shows that </>rlp depends on the frictional properties of 
the beads, as predicted in 

Conclusions. We have shown that the configuration 
space of a granular medium can, by using flow pulses, be 
explored in a history independent way, which is essential 
for the statistical theory [l| to be applicable. Fluctua- 
tions of the volume fraction are Gaussian with a parabolic 
minimum, which corresponds to a maximum in the num- 
ber of statistically independent regions. The minimum 
arises as a consequence of competing mechanisms; the 
location of the minimum depends on the frictional prop- 
erties of the beads. The compactivity can be determined 
from the volume fraction fluctuations; hence compactiv- 
ity is a well-defined parameter representative of the mate- 
rial. This opens the door for new experiments that would, 
for example, investigate the use of compactivity as a con- 
trol parameter for segregation in binary mixtures 0, Il5| . 
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